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Abstract: Double hard scattering in proton-proton collisions is described in terms of 
double parton distributions. We derive bounds on these distributions that follow from their 
interpretation as probability densities, taking into account all possible spin correlations 
between two partons in an unpolarized proton. These bounds constrain the size of the 
polarized distributions and can for instance be used to set upper limits on the effects of 
spin correlations in double hard scattering. We investigate the stability of the bounds 
under leading-order DGLAP evolution to higher scales. 
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1 Introduction 

In a time when the dynamics of the strong interaction in hadron-hadron collisions is mov- 
ing towards the domain of precision physics, there are still aspects that are under poor 
theoretical and experimental control. One of these aspects is double parton scattering, 
where two partons from each proton have a hard interaction in a single proton-proton 
collision. Correlations between the two hard interactions have been the subject of several 
recent studies [1-8]. The relevance of spin correlations in double parton scattering was 
pointed out long ago [9, 10] and recently followed up in [11, 12]. The studies in [13] and 
[14] have shown that spin correlations in the production of two vector bosons by double 
hard scattering have observable effects both on the interaction rate and on kinematic dis- 
tributions. Spin correlations between the two partons are quantified by polarized double 
parton distributions (DPDs), which describe for instance the difference of the probability 
densities for finding two quarks with equal or with opposite helicities. It was argued in [12] 
that such correlations need not be small, and a recent study in the MIT bag model [15] 
indeed found large spin correlations between quarks in the valence region. However, our 
knowledge of polarized DPDs is still poor at best, and any information about them is of 
value. 
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In the present work, we derive model independent constraints on DPDs that follow 
from their interpretation as probability densities for finding two partons in a specified po- 
larization state. Similar positivity bounds have been derived for single-parton distributions 
in the form of the Soffer bound [16] and of inequalities for transverse-momentum dependent 
distributions [17] and generalized parton distributions [18]. 

The structure of this paper is as follows. In the next section we set the stage by 
introducing the DPDs for different polarizations and parton species. In section 3 we derive 
the spin density matrices for two partons inside an unpolarized proton, and in section 4 
we use these matrices to derive bounds on polarized DPDs. In section 5 we show that 
the homogeneous leading-order evolution equations preserve these bounds when going to 
higher scales. We conclude in section 6 and give some technical details in two appendices. 



2 Double parton distributions 

Double parton distributions for quarks and antiquarks have been extensively studied in 
[12], and we only review the properties important for our purpose. Since we will need a 
probability interpretation, we restrict ourselves to distributions that are integrated over 
the transverse parton momenta and that have a trivial color structure. In the parlance of 
[12] these are collinear color-singlet distributions. 

Collinear DPDs depend on the longitudinal momentum fractions x\ and x 2 of the two 
partons and on the transverse distance y between them. For two partons a\ and a% in an 
unpolarized right-moving proton we write 

F aia2 (x 1 ,x 2 ,y) = 2p + (x lP + )-^ (x 2 p + )-" 2 J^-^-dy- e ^T+^)P + 

x ip\O a2 (0,z 2 )Oa 1 (y,z 1 )\p) , (2.1) 

where ni = 1 if parton number i is a gluon and = otherwise. We use light-cone 
coordinates v ± = (v° db v 3 )/\^2 and the transverse component v = {v l ,v 2 ) for any four- 
vector v. The operators for quarks read 

O ai {y,z i ) = q i (y-\zi)T a .q i (y + \zi) (2.2) 

*i=V -0, Z;=0 

with projections 

r g = h + > r Ag = i 7 + 75 , Y j Sq = \ia^ l5 (j = 1,2) (2.3) 

onto unpolarized quarks (q), longitudinally polarized quarks (Aq) and transversely polar- 
ized quarks (8q). The field with argument y + \zi in O ai (y, zi) is associated with a quark 
in the amplitude of a double scattering process and the field with argument y — \zi with 
a quark in the complex conjugate amplitude. The operators for gluons are 

O ai {y,z i ) = Wj:G^'{y-\z i )G^{y+\z i ) + (2.4) 

«i =!/ -0, z i= 
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with projections 

Uf = 6 jj ', n%' g = ie jj ', [ng'P" = T jj '' kk ' (2.5) 

onto unpolarized gluons (g), longitudinally polarized gluons (Ag) and linearly polarized 
gluons (5g). The tensor 

T jj',kk' = i ( § jk 6 j>k> + pk'p'k _ 6 jj' § kk>^ ^ 

satisfies T oi ,kk T kk,ll _ T jj ,ll an j j g S y mm etric and traceless in each of the index pairs 
(jf) and (kk'). Note that for gluons 5g denotes linear polarization, i.e. the interference 
between gluons whose helicities differ by two units in the scattering amplitude and its 
conjugate, while for quarks 5q symbolizes transverse polarization, where the interference 
is between quarks with a helicity difference of one unit. Since we limit ourselves to color- 
singlet distributions, a sum over the color indices of the quark fields in (2.2) and the gluon 
fields in (2.4) is implied. We do not write out the Wilson lines that make the operators 
gauge invariant. 

The different spin projections lead to a large number of DPDs. For collinear color- 
singlet distributions, several polarization combinations are zero due to time reversal and 
parity invariance. This concerns the DPDs with one longitudinally polarized and one 
unpolarized parton, as well as those with one longitudinally polarized parton and one 
transversely polarized (anti)quark or linearly polarized gluon. A decomposition of the 
nonzero distributions for two quarks in terms of real-valued scalar functions has already 
been given in [12]: 

F qq (x 1 ,X 2 ,y) = f qq (x 1 ,X 2 ,y) , 
FAqAq{xi,X 2 ,y) = fAqAq{xi,X 2 ,y) , 

F 3 qS q{xi,x 2 ,y) = y 1 Mf qSq (x 1 ,x 2 ,y) , 
F s qq (xi,x 2 ,y) = y :, Mfs qq {xi,x 2 ,y) , 

F^ 5q ( Xl ,x 2 ,y) = S^'f SgSq (x 1 ,x 2 ,y) + 2^'™ M 2 f SqSq (x u x 2 , y) , (2.7) 

where M is the proton mass, yi = e JJ ' y J ' and y = \p\p ■ We use a shorthand notation where 
vectors y or y appearing as an index of r denote contraction, i.e. t^' ,vv = T^' ,kk ' y k y k 
etc. Decompositions analogous to (2.7) hold for quark-antiquark distributions and for the 
distributions of two antiquarks. 

Since quarks and gluons mix under evolution, we also need to consider DPDs involving 
gluons. We define 

F qg {x 1 ,x 2 ,y) = f qg (x 1 ,x 2 ,y) , 

FAqAg{xi,X 2 ,y) = fAqAg(x 1 ,X 2 ,y) , 

F 3 J g ( Xl ,x 2 ,y) = T^M 2 f qig {x^x 2 ,y) , 
F 5 qg {xi,x 2 ,y) = y J Mfs qg (x 1 ,x 2 ,y) , 
Fl q fg\xi,x 2 ,y) = - TV> kW Mf Sq s g {x x ,x 2 ,y) 

_ ( y j T kk',yy + y j T kk>yy} M 3 /| g5fl (^l, X 2 , y) (2.8) 
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for quark-gluon distributions, with analogous expressions for gluon-quark distributions and 
distributions where the quark is replaced by an antiquark. For two-gluon distributions we 
write 



We remark that, although linear gluon polarization is described by a tensor with two 
indices, the restriction that this tensor is symmetric and traceless gives rise to the same 
number of distributions as for transverse quark polarization, which is described by a vector. 
The prefactors in (2.8) and (2.9) have been chosen such that we will obtain a simple 
correspondence between quark and gluon distributions in the spin density matrices to be 
derived in the next section. 

Note that DPDs involving gluons are not only relevant in the context of evolution but 
also enter directly in important double scattering processes such as the production of jets. 
Their properties are hence of considerable practical interest. 

In complete analogy to the case of collinear single-parton distributions, the DPDs we 
have introduced can be interpreted as probability densities for finding two partons inside an 
unpolarized proton, with a relative transverse distance y and with longitudinal momentum 
fractions x\ and X2- This becomes evident from their appearance in the cross section 
formulae for double parton scattering [12]. It can also be seen from a representation in 
terms of parton creation and annihilation operators or from a representation in terms of 
the light-cone wave functions of the proton, which are straightforward extensions of the 
corresponding representations for single-parton distributions (given for instance in sections 



As in the case of single-parton densities, this interpretation does however not strictly 
hold in QCD, because the distributions are defined with subtractions from the ultraviolet 
region of parton momenta. The subtraction terms can in principle invalidate the positivity 
of the distributions. Nevertheless, it is useful to explore the consequences of the probability 
interpretation as a guide for developing physically intuitive models of the distributions. 
This holds in particular if one works in leading order of a s , where the connection between 
parton distributions and physical cross sections (which must of course be positive semi- 
definite) is most direct. 

3 Two-parton spin density matrices 

The polarization state of two partons in an unpolarized proton is described by a spin 
density matrix that can be written in terms of the DPDs we introduced in the previous 





f gg (x 1 ,x 2 ,y) , 

fAgAg(xx,X2,y) , 



jj',yy--kk', 



)M 4 f SgSg (x 1 ,x 2 ,y). (2.9) 



3.4 and 3.11 of [19]). 
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section. We start by trading the projection operators (2.3) and (2.5) for operators that 
project onto quarks or gluons of definite helicity. We can then easily write down the spin 
density matrix for two partons in the helicity basis. 

The projection operators T\>\ for quarks, where A (A') refers to the quark helicity in 
the amplitude (conjugate amplitude), are given by 



r++ 



7 



(1+75) 
(1-75) 



r 
r 



ia 



+i 



(1 + 75) 



r 1 



iV 2 
11 Sq 



(3.1) 



Here we use the phase conventions for spin-half fields specified in [20]. The projection 
operators Ily* for gluons, where A and j (A' and j') refer to the amplitude (conjugate 
amplitude), can be constructed from the polarization vectors 



and read 



V2 



:m). 



it" 
1L ++ 



U 33_ 

u jj ' 



i 

7^ 



1 ( rrJ'j" -i- TT-"' 
z ( Txji' _ t\H' 

2 \ro A 9 

rill 33' 



i,-0 



(3.2) 



(3.3) 



We can now organize the distributions in matrices where the columns (rows) correspond to 

helicity states ++, — h, H — , of the two partons in the amplitude (conjugate amplitude). 

The spin density matrix for two quarks reads 



' fqq + fAqAq 

ie-^yyMfs, 



-ie^yyMf 5qq -ie^yMf qSq 2e 2i ^y 2 M 2 fl qSq \ 



qq 

ie-^yyMf q5q 



fqq — fAqAq 
2f6qSq 



2f6q5q 
fqq — fAqAq 



-ie^yyMf qS q 
-ie^yyMfs, 



qq 



(3.4) 



2e~ 2 ^y y 2 M 2 f 5qSq ie-^yyMf q Sq ie-^yyMfsqq f qq + f Aq A q j 



where the angle ip y describes the orientation of the vector y = y (cos <p y , sin (p y ) in the 
transverse plane. With the caveat spelled out at the end of the previous section, the diago- 
nal matrix elements can be interpreted as the probability densities for finding two partons 
in definite helicity states inside an unpolarized proton. Specifically, f qq + fAqAq is the 
probability density for finding two quarks with positive helicities, which in an unpolarized 
proton is equal to the probability density for finding two quarks with negative helicities. 
The probability density for finding two quarks with opposite helicities is f qq — fAqAq- The 
off-diagonal elements of p describe helicity interference, with f$ q g q in the right upper corner 
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corresponding for instance to the case where both quarks have negative helicity in the am- 
plitude and positive helicity in the conjugate amplitude. This leads to a helicity difference 
between the amplitude and its conjugate, which is balanced by two units of orbital angular 
momentum indicated by an exponential e 2%LPy and an associated factor y 2 . By contrast, 
fsqSq describes the case when the helicity difference is +1 for one quark and —1 for the 
other, so that the overall helicity is balanced. 

Turning now to gluons, we have a spin density matrix 



/ fqg + fAq 



IqgT j AqAg -%e 1 ^ yMf Sqg 

ie-^yMfsqg f qg - fA q A g 
. e -2i Vyy 2 M 2 fqgg 2ie-^yyMf SqSg 



-e 2il ?yy 2 M 2 f qS g -2ie 3i ^y 3 M 3 /, 



SqSg 



-2ie^yyMfs qS g 

fqg ~ fAqAg 



^2ze- 3i ^y 3 M 3 / ( 



t 

SqSg 



-2iip l 



y 2 M 2 f q5g ie-^vyMfs, 



</.'/ 



e 2 ^yy 2 M 2 f q5g 
-ie^yyMhqg 

fqg + fAqAg 



(3.5) 



for quark-gluon distributions and an analogous matrix for gluon-quark distributions. For 
two-gluon distributions we find 



' fgg + /a S A 9 

-e- 2 ^yy 2 M 2 f 5gg 
-e- 2i ^y 2 M 2 f gSg 
\—2e~ 4i<p vy 4 M 4 f ( 



-e 2l ^y 2 M 2 f Sgg 

fgg ~ fAgAg 
2f&gSg 



D 2iip. 



•y 2 M 2 f gSg -2e 4 ^y 4 M 4 /L g \ 



2f SgS g 
fgg ~ fAgAg 



,2i<p yy 2 M 2 fgSg 

z 2iLp yy 2 M 2 f Sg g 



-e- 2i ^y 2 M 2 fg S g -e- 2 ^yy 2 M 2 f 5g g f gg + f AgAg ) 



(3.6) 



The matrices for distributions where quarks are replaced by antiquarks are analogous to 
(3.4) and (3.5). We see that the parameterization of DPDs in the previous section gives 
simple expressions for the spin density matrices and similar structures for all types of 
partons. 

The difference in spin between quarks and gluons causes the different dependence on 
the azimuthal angle ip y in (3.4), (3.5) and (3.6). A mismatch of n units between the sum of 
parton helicities in the amplitude and its conjugate goes along with an exponential e ±niip y 
and an associated factor y n . 



4 Positivity bounds 

We now show how the probability interpretation of DPDs constrains the size of the polarized 
distributions. Since the probability density for finding two partons in a general polarization 
state is positive semi-definite, we have 

v l' 1 \' 2 P(x' 1 \' 2 ){x 1 \ 2 ) v x 1 \ 2 ^ 

X' 1 X' 2 X 1 X 2 

with arbitrary complex coefficients v\ 1 \ 2 normalized as X^AiA 2 l t, AiA 2 | 2 = 1- The helicity 
matrices are therefore positive semi-definite. The same property has been derived for the 
spin density matrices associated with transverse-momentum dependent distributions [17] 
or generalized parton distributions [18]. 
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To simplify the algebra, we first cast all helicity matrices into a common form that is 
independent of the angle ip y . This is achieved by unitary transformations, multiplying by 
a matrix U from the right and by W from the left. The transformation matrices for the 
parton combinations in (3.4) to (3.6) are 



U qq = diag(-e 2i< ^, -ie i<fiy , - 



U qg = diag(ze 



3iip y p 2iifiy 



Itfy 



u, 



99 



diag(e 4 ^,-e 



2iip y 



!)■ 



(4.2) 



After these transformations and their analog for gluon-quark distributions, the spin density 
matrices can be written as 



/ fab + fAaAb hs a b 



h 



aSb 



-Han ^ 



V 



hsab 
h a Sb 

-9/7* 

zn Sa8b 



fab — fAaAb ^hgaSb h a Sb 
ZhsaSb fab ~ fAaAb hs a b 
h a 8b hsab fab + fAaAb ) 



(4.3) 



with the following identification of distributions for different parton combinations: 

fab = fqq 5 fqg j fgq 5 fgg > 
fAaAb = fAqAq , f AqAg , f AgAq , fAgAg , 

hsab = yMfs qq , yMf Sqg , y 2 M 2 f Sgq , y 2 M 2 fs gg , 



Ksb = yMfgSg , y 2 M 2 f qog , yMfgSq , y l M l f g0 g 



2 7\ /r2 



hsaSb = fSqSq , yM fsqSg > yMfSgSq , fSgSg , 



n SaSb 



y 2 M 2 f SqSq , y 3 M 3 fl qSg , y 3 M 3 f Sg5q , y*M*fl gSg 



(4.4) 



Analogous expressions hold if quarks are replaced by antiquarks. Positivity 1 of the diagonal 
elements of p yields the trivial bounds 



fab > /AaAb 



(4.5) 



The principal minors of the two-dimensional sub-spaces must be positive semi-definite as 
well, which gives upper bounds on the distributions for one or two transversely or linearly 
polarized partons: 

fab + fAaAb > 2 I h l SaSb I , 



fab — fAaAb > 2 | hfiaSb \ , 
fab > (fab + fAaAb)(fab ~ fAaAb) > hf 



Sab 1 



fab - (fab + fAaAb)(fab ~ fAaAb) > h aSb . 



(4.6) 



The principal minors of dimension three, as well as det(p) provide further bounds, which 
are rather cumbersome and will not be given here. The strongest bounds can be obtained 



For ease of language we use "positivity" in the sense of "positive semi-definite" here and in the following. 
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from the positivity of the eigenvalues of p, which is a sufficient and necessary condition for 
the positivity of p. Calculating the eigenvalues we obtain 

fab + hgaSb - h\ a5b ± J (h Sab + h aSb ) 2 + (/AaAfe ~ hs a Sb ~ Ma<5b) 2 - ' 

fab ~ hsaSb + hgaSb ± \j {hSab ~ Ksb) 2 + (/AaAfe + h SaSb + h^) 2 > . (4.7) 

These inequalities set upper limits on the size of spin correlations between two partons in 
an unpolarized proton. They can be used either to construct double parton distributions 
or to put limits on polarization effects in double hard scattering processes. 

We note that positive semidefinite combinations of DPDs were discussed already in 
the pioneering studies [9, 10]. Distributions that involve a helicity mismatch between the 
amplitude and its conjugate (see section 3) were however not considered in that work. 
The derivation in [9, 10] thus corresponds to our results (4.5) and (4.6) if all distributions 
multiplied with a power of y in (4.4) are set to zero. 

5 Stability under evolution 

The ultraviolet subtractions mentioned at the end of section 2 induce a scale dependence, 
which for collinear single-parton distributions is described by the DGLAP evolution equa- 
tions. While the subtractions themselves may invalidate positivity of the distributions and 
thus their density interpretation, the evolution equations can be interpreted in a proba- 
bilistic manner provided that one takes the leading-order approximation of the evolution 
kernels [23, 24]. Specifically, one finds that if parton distributions are positive semi-definite 
at a certain scale, this property is preserved by leading-order evolution to higher scales. 
This also holds for the Soffer inequality, which expresses positivity in the sector of trans- 
verse quark polarization [25, 26]. For evolution at next-to-leading order in a s the situation 
is less clear-cut and a discussion of positivity depends in particular on the scheme in which 
the distributions are defined. We refer to [27] and [28, 29] for a discussion of the situation 
for longitudinal and transverse parton polarization, respectively. 

Returning to double parton distributions, we now show that the bounds derived in the 
previous section are stable under leading-order evolution to higher scales. The strategy 
for the derivation is as follows: we first introduce linear combinations of double parton 
distributions whose positivity is necessary and sufficient for the positivity of the spin density 
matrices and then show that these linear combinations remain positive semi-definite under 
evolution. The positivity of the spin density matrices then guarantees the stability of the 
positivity bounds. 

5.1 Evolution of double parton distributions 

To begin with, let us specify the evolution of collinear DPDs in the color-singlet sector. 
We use the homogeneous evolution equations in the transverse position representation, 
see e.g. equation (5.93) in [12]. These equations apply at nonzero y if F aia2 (xi,X2,y) is 
defined via (2.1) with the operators O ai (y,zi) and O a2 (0, Z2) renormalized by standard 
MS subtraction. The inhomogeneous term for the splitting of one parton into two that 
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has been previously considered in the literature [30-34] does not appear in this case. As 
discussed in [11, 12], a consistent formulation of factorization for double parton scattering 
does not yet exist, so that it remains unclear how DPDs should best be defined (and how 
they evolve). For simplicity we will limit our present investigation to the homogeneous 
evolution equations. 

It is useful for our purpose to take different renormalization scales fi\ and fi2 for the 
two partons, corresponding to separate ultraviolet renormalization of O ai and O a , 2 in (2.1). 
The evolution equation for the unpolarized double quark distributions in the first scale 
then reads 



df qq (x 1 ,x 2 ,y,Hi,fJ>2) 
dri 



Pqq®l fqq + Pqg®l fgq, (5.1) 



where 



pl—X2 / x±\ 

J XI u l V"l/ 

is a convolution in the first argument of the DPDs with the leading-order splitting functions 
P a b known from DGLAP evolution of single-parton distributions. We note that the leading- 
order splitting functions are the same for quarks and antiquarks, i.e. one has P qq — Pqqi 
P qg = P qg , P gq = P gq and analogous relations for polarized partons. In appendix A we 
give the explicit evolution equations for all polarized DPDs and list the associated splitting 
functions. 

The evolution variable in (5.1) is taken as 

T1 = J — > (5 - 3) 

where the lower limit of integration is irrelevant in the derivative df /dr\. The use of t\ 
is just a matter of convenience as it removes the running coupling from the leading-order 
splitting functions. 

The analog of (5.1) for the scale associated with the second parton is 



dfqq(x 1 ,X2,y;(J,l,(J-2) 



Pqq ®2 fqq + Pqg ®2 fqg ■ (5.4) 



<9r 2 

The evolution equation for equal scales, i.e. for f qg (xi, X2, y; /U, /x), is readily obtained by 
adding the right-hand sides of (5.1) and (5.4). We will show that positivity is preserved 
for separate evolution in fix. The same then obviously holds for evolution in [12 and hence 
for the evolution in a single common scale [i\ = [12- 

5.2 Linear combinations of DPDs 

A key ingredient in our argument is to form suitable linear combinations of double parton 
distributions, which we now introduce. Positivity of the spin density matrix p means that 
v'pv > for any complex vector v, as we spelled out in (4.1). Parameterizing the vector 
as 

v T = (ai + i bi , a 2 + ib 2 , a 3 + ib 3 , 04 + ib±) (5.5) 
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with real numbers aj, 6, and performing the multiplication with the matrix in (4.3) gives 

Q~a~b = c ^fab + C 2 h a Sb + C 3 fAaAb + C^h oab + C 5 h Sa 8b + C^^SaSb > i ( 5 - 6 ) 

where = 4t^/9 u and the coefficients q are given by 

00000000 

ci = a 1 + b 1 + a 2 + b 2 + a 3 + b 3 + a A + b A , c 2 = 2(aia 3 + 6163 + a 2 a 4 + &2&4) , 
c 3 = + 6f — al — b\ — a| — b\ + a| + b\ , c 4 = 2(aia 2 + b\b 2 + a 3 a 4 + 6364) , 
C5 = 4(a 2 a 3 + 6 2 6 3 ) , c 6 = -4(oio 4 + 6164) . (5.7) 

We will prove the stability of the positivity bounds by showing that for arbitrary values of 
en and bi the inequality (5.6) is stable under evolution to higher scales. It will be convenient 
to consider further linear combinations of double parton distributions. Changing signs of 
the parameters a± — > — a\, b\ — > —bi, a 3 — > — a 3 , 6 3 — > —63 we get 

Qab = c ^fab + C 2 h aS b + C 3 f AaAb ~ C4 hs a b ~ C5 hs a Sb ~ CQ h\ a Sb > • (5.8) 

Adding (5.6) and (5.8) gives the simpler inequality 

B ab = c ^fab + c 2 h aSb + c 3 f AaAb > , (5.9) 
and interchanging indices (lf>2 and 3 <B- 4) in the elements of v gives 

B ab = c ^ob + c 2 h aSb - c 3 f AaAb > . (5.10) 

If (5.6) holds at a given scale for arbitrary values of o« and bi, then (5.8) to (5.10) hold at 
that scale as well. 

We will see that the evolution equations in the scale \x\ can be formulated in terms of 
Qab' Qab an< ^ B ab a l° ne - 2 This becomes plausible if we note that these three functions are 
linear combinations of {cif ab + c 2 h a5b ), f Aa Ab and (c4hsab + C5h Sa s b + c 6 h t SaSb ) and that for 
evolution in fi\ only the polarization of the first parton is relevant but not the polarization 
of the second parton. The linear combinations Q^ b may be regarded as generalizations of 
the distributions Q± = \{q + q) ±5q introduced in [26], where it was shown that the Soffer 
bound for the quark transversity distribution 5q is stable under leading-order evolution to 
higher scales. 

5.3 Evolution of the linear combinations 

We now show that the distributions Q^ b and B^ h remain positive semi-definite under 
leading-order evolution to higher scales. This implies the positivity of the spin density 
matrices and thereby the validity of the bounds derived in section 4. 
The evolution equations for the distributions Q^ b are 

VP q - q sp q + q Mq-J \p+pJ 






V 99 qq 



2 The combination B^ b = (Q^ b + Q ab )/2 is not independent and just used as an abbreviation. 
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for a quark as first parton and 




5P + 5P~ 

u 99 99 

5P~ 5P + 

VI gg VI gg 



99 99 
p- p- 
K 99 99 




p+ p- 

ga ga 
p+ p- 
ga ga j 



when the first parton is a gluon. The evolution equations for B , read 




for a quark and 



p+ p-< 
n n 
p- p+ 
K qq qq. 




P + P -' 

99 19 
p- p+ 

19 qg / 



gb 
gbi 




(5.12) 



(5.13) 




P + p- 
gg gg 
p- p+ 
v gg gg, 



gb 
gb/ 



+ E 



a=q,q 



p+ p- 
ga ± ga 

Pga Pga. 



B 



ab 
ab i 



(5.14) 



for a gluon. The evolution equations have the same form for antiquarks, i.e. (5.11) and 
(5.13) remain valid if we replace q — > q everywhere (except in the label b for the second 
parton, which always remains fixed when we consider evolution in fi\). 

The splitting functions appearing in the above equations are defined as 



3± 



1 



{Pab ± PAaAb) , 



1 



SP a b = n {P&aAb ± PSaSb) 



(5.15) 



for all parton indices a and b. We remark that the kernels P^ b {P~ h ) correspond to the case 
where the parton helicity is conserved (flipped). The only splitting functions that receive 
contributions from virtual graphs and hence contain a plus-prescription or an explicit 5 
function are 



P 



Cf 
2 



5P + = ^ 
ii 2 



P, 



gg 



2N C 
5P; g = 2N C 
with N c = 3,C F = 4/3 and 



2(1 + z 2 ] 
(1-*)+ 
(1 + ^) 2 

(1-*)+ 

z 

z 



- + 3(5(1 - z 
+ 3S(l-z) 



+ 



(l-*)(l + z) s 



(1 -z) 



+ 1 



+ 



2.; 



+ — 6{l-z), 



11 2 

Pa = y Nc ~ 3 n /' 



(5.16) 



(5.17) 



where rif is the number of active quark flavors. They are all positive for < z < 1 but 
have negative contributions at z = 1 that arise from the plus-prescription, whose form 
is recalled in (B.3). In appendix B we show explicitly that the virtual contribution to 
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evolution cannot change the sign of the distributions, which has previously been argued 
to be the case based on the probabilistic interpretation of leading-order evolution and its 
relation to the Boltzmann equation [23, 24, 26]. The reason for this property is that the 
virtual contribution to the evolution of a function is proportional to the function itself. 
We can then conclude that the diagonal terms in the evolution equations (5.11) to (5.14) 
preserve positivity. The off-diagonal kernels 

Pa~ q = 0, P~ g = N c {l-zf/z, 

SP- =C F (1- z)/2 , 5P~ = 2N C (1 - z) (5.18) 

and 

^ = ^/2, P+ = C F /z, 

P- = (1 - zf/2 , Pgq = C F (1- zf/z . (5.19) 

are all positive or zero for < z < 1 and regular at z = 1. Therefore they only reinforce 
positivity. In summary, if we have positive semi-definite initial conditions for all functions 
Q^ b and B^ b at some scale, then evolution to higher scales preserves this property. A more 
explicit derivation is given in appendix B. 



6 Conclusions 

We have derived spin density matrices for double parton distributions of quarks, anti- 
quarks and gluons. These matrices reveal the full polarization structure of two partons in 
an unpolarized proton and show the correspondence between the different polarized double 
parton distributions and parton helicities. The probabilistic interpretation of the double 
parton distribution for an arbitrary polarization state of the two partons gives upper limits 
on the size of spin correlations. These positivity bounds can be useful for modeling the 
otherwise poorly constrained double parton distributions and for deriving upper limits on 
spin effects in double hard scattering processes. We have shown that the homogeneous 
leading-order evolution equations preserve the bounds when going from lower to higher 
scales. 



A Evolution equations and splitting functions 

For completeness we give here the leading-order evolution equations for the first parton in 
the double parton distributions. When the first parton is a quark, we have 

Pqq ®1 fqb + Pqg ®l fgb i 



dfAqAb 



= Pqq ®1 fqSb + Pqg ®1 fgSb , 



PAqAq <8>1 fAqAb + P&qAg ®l fAgAb ■ 



dn 

dfsqb _ D „ , QfsqSb _ _ . dfsqSb 



= PSqSq ®1 fsqb , = P&q5q ®1 fSqSb , = PSqSq ®1 fsqSb (^--1) 
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for b = q,q,g. The arguments of the distributions are as in (5.1) and (5.2). Analogous 
equations hold if the first parton is an antiquark. For gluons we have 



dfgb 
QfgSb 



Pgg ®1 fgb + ^2 Pga <8>1 fab , 
a=q,q 

Pgg ®1 fgSb + ^ P 9 a ® l ^ a5 



a=q,q 



"^f Ajl = PAgAg ®1 /AgAb + ^ -PA s Aa <8>1 /AaAfe , 



3t\ 

dfsgb 



8t\ 



Pi 



1 fSgb i 



9ri 



-P5 9 5 9 ®1 f&gSb , 



df, 



8t\ 



P S g5g ®1 fsgSb ■ (A.2) 



The leading-order splitting functions have been derived in [21, 22]. They are given by 

l + z 2 



P qq (z) = C F 

PAqAq{z) = Pqq(z) , 
P 5 q8q(z) = Pqq(z) - C F (l - Z) 

for quark-quark transitions and by 

(1-zKl + z*) 



(A.3) 



Pgg(z) = 2N C 



PAgAg(z) = Pgg(z) - 2N C 
P S g5g(z) = Pgg(z) - 2N C 



1 - z)+ + 

(1-zf 



+ — 5(l-z) 



(l-z)(l + z 2 ) 



(A.4) 



for gluons. The splitting functions that mix quarks and gluons read 



P 



19 



P AqAg 



z 2 + (l-z) 



:i-z) 



Pgq - C F 



l + (l-z) 



PAgAg — Cf 



\-{\-z\ 



(A.5) 



As already mentioned below (5.2), the splitting functions are identical for quarks and 
antiquarks, i.e. (A.3) and (A.5) remain valid if we replace q — > q. At leading order in a s 
there are no direct transitions between quarks and antiquarks. 



B Elements of a stability proof 

In this appendix we show in more detail that the evolution equations in section 5.3 preserve 
positivity, taking particular care of the negative terms in the splitting functions that arise 
from virtual graphs and are implicit in the plus-prescription. We first consider the evolution 
of a single distribution and then extend the argument to the full coupled system of evolution 
equations. 
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We examine a function evolving as 



a 



v du 



-/(x,r)= / -P(-)f( u ,r 



u 



u 



with < x < v < 1 and separate the splitting function as 

P(z) = jp^r- + P r (z) + P S 6(1 - z) , 
v 1 - z )+ 



(B.l) 



(B.2) 



where P s {z) and P r (z) are positive semi-definite for < z < 1 and regular at z = 1. The 
constant may be positive, negative or zero. The plus-prescription is defined as usual by 



[s(z)} + = s(z) - 5(1 -z) dz' s(z') 



(B.3) 



where it is understood that the non-integrable singularity in the last term cancels when 
(B.3) is integrated over with a smooth test function. The plus-prescription part of the 
convolution in (B.l) can be written as 



du P s (x I u) 
u (1 — x/u)+ 



f(n,r) 



" P s (x/u) 

du j (u, T) + 

X+E 



u — X 







du^^f(x,r) + 0(e) 
u — X 



(B.4) 



where for the error estimate we have assumed that f(u, r) is differentiable at u = x. 
Defining 



9e(x,T-J) 

h e (x) 



du 



P s (x/u) P r (x/u) 



Ps + Ps(l) 



u — X u 
" x ~ e du 



f(u,r) 







x — u 



we can approximate the evolution of / by 







f(x, t) = g e (x, r; /) - h e (x) f(x, r) 



(B.5) 



(B.6) 



with an error that becomes arbitrarily small for e — > 0. In a more formal proof, one would 
replace / with f t in (B.6) and show that lim f € is a solution of (B.l) . We now rewrite 
(B.6) as 

d r 



e->0 



dr 



^f(x,r) =e Th ^g € (x,r;f) 



(B.7) 



Since g 6 is the convolution of f(x,r) with a positive semi-definite integral kernel, the r.h.s. 
of this equation is positive semi-definite as long as f(x,r) is. With initial conditions 
f(x,To) > for all x at a starting scale To, the function e Thc ^ x ' ) f(x, r) can therefore not 
decrease as r increases, so that /(x,r) stays positive semi-definite for all r > tq. We note 
that the sign of h e (x) and thus of the constant Pg is irrelevant for this argument. 
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We now consider the coupled system of evolution equations given by (5.11) to (5.14). 
Using a vector notation f l (x, r) for the 8ny + 4 functions Q^ b , Q~ b , B^ b , B~ b with a = q,q,g 
(and b fixed), we can cast their evolution into the form 

Y f(x, t) = gl(x, r; f) - hi(x) f(x, r) + ^ & Q f j (% r) (B.8) 

T . + . J x U 

with i = 1, . . . , 8nf + 4. Here g\ and h\ are defined as in (B.5) with regular and positive 
semi-definite functions Pl{z) and P^iz). The mixing kernels P l;, (z) in (B.8) are regular 
and positive semi-definite as well. Rewriting the evolution as 



e Th ^f(x,T) 



Ot 



e rh E (x) 



(B.9) 



we see that if one has initial conditions f 3 {x,To) > for all j then all functions f J (x,T) 
remain positive semi-definite for r > tq. 
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